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The bifurcation period in low-mass X-ray binaries is the initial orbital pe- 
riod which separates the formation of converging systems (which evolve with 
^ decreasing orbital periods until the donor becomes degenerate) from the diverg 

O 

I atically the bifurcation periods of binary systems with a IAMq neutron star and 
a 0.5 — 2Mq donor star, taking into account different kinds of magnetic brak- 

^ I ing and mass loss mechanisms. Our results show that the saturated magnetic 

Q\ \ braking can considerably decrease the values of bifurcation period compared to 

^ I the traditional magnetic braking, while the influence of mass loss mechanisms on 

CN I bifurcation periods is quite weak. We also develop a semi-analytical method to 

O I compute the bifurcation period, the result of which agrees well with the numerical 

00 ' method in the leading order. 
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ing systems (which evolve with increasing orbital periods until the donor star 
loses its envelope and a wide detached binary is formed). We calculate system- 



Subject headings: binaries: close - stars: evolution - X-rays: binaries 



1. Introduction 

One interesting and important topic in the secular evolution of low-mass X-ray binaries 
(LMXBs) is the so-called "bifurcation period" Pbif, the initial binary orbital period which 
separates the formation of converging systems (which evolve with decreasing orbital periods 
until the donor becomes degenerate) from the diverging systems (which evolve with increas- 



ing orbit al periods unti 



formed) (ITutukov et al 



the donor star loses its envelope and a wide detached binary is 



19851). The first system atic investigations on the bifurcation period 



were done by iPylyser fc Savonijd (119881 . Il989l ) . Neglecting mass loss fr om the binary sys 



tem and assuming angular momentum loss due to magnetic braking (MB; IVerbunt fc Zwaan 
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198ll ) and gravitational radiation (GR; iLandau &: Lifshitall975l ). these authors found that 
the bifurcation period is in the range Phi f ^ 0-4 — 0.7 d ay for LMXBs, and strongly de- 



pends on magnetic braking efficiency. lErgma et al.l (119981 ) included mass loss from the bi- 
nary system and re-calculated the bifurcation period for two mass configurations (Mi/Mq, 
M2/M0) = (1.4,1) and (1.4,1.5) and two chemical compositions {Z = 0.003, 0.03). They 
pointed out that the mass loss from the binary system also plays an impotent role besides 
magnetic braking in determining the value of Pbif, while the chemical composition could 
only cause small change in Pbif- Their bifurcation periods are Pbif ~ 0.85 — 1.05 day under 
conservative mass t ransfer, and 1.6 — 1.7 time s larger if moderate non- conservative mass 
transfer is assumed. iPodsiadlowski et al.l (120021 ) found a bifurcation period around 18 hr for 
a I.4M0 NS and a IMq companion star, where they defined the bifurcation period as the 
orbital period when the Roche lobe overfiow just began, instead of the initial orbital period. 



van der Sluys et al.l (l2005al Jbl) also investigated the bifurcation period in LMXBs fo- 
cusing the formation of ultra-compact X-ray binaries (UCXBs), and specified the bifurca- 
tion period as "the longest initial period that leads to UCXBs within a Hubble time (13.7 
Gyr)". UCXBs are bright X-ray sources with very short orbital periods (P < 1 h). The 
donor has to be a compact source like a white dwarf or a compact core of an evolved giant 
star to fit in the small Roche lobe size. Such sources may be formed through dynami- 



cal processes including stellar collisions and common envelop evolution (j Clark et al.lll975 



Rasio et al.l I2OOOI : iLombardi et al.l l2006l ) . An alternative scenario for the formation of such 
sources is through stable mass transfer in X-ray binaries with a low- or intermediate-mass 
donor star, which may explain the negative der ivative of the 11-min source in NGC 6624 
( Ivan der Klis et al.lll993l : IChou fc Grindlayll200ll ). It has been found th at systems with ini- 



tial orbital period just below the bifurcation period may form UCXBs (iNelson et al.l 



Tutukov et al. 
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19871: 



Pvlvser &: Savoniielll988l : IPodsiadlowski et al.ll2002l : Ivan der Sluys et al. 



1986 



2005al ). IPodsiadlowski et al.l (120021 ) showed that the closer the initial orbital period to the 
bifurcation period from below, the smaller the minimum orbital period will be achieved. 
So the value of bifurcation period is crucial to understanding the formation of UCXBs 
( Ivan der Sluys et al.ll2005bl ). 



In this paper we make a systematic investigation on the bifurcation period for binary 
systems containing an NS with a main-sequence (MS) companion of mass from 0.5Mq to 
2Mq. This work was motivated by recent progress in studies on mass and angular mo- 
mentum loss mechanisms in LMXB ev oluti on. In previous works the MB law originally 
postulated by IVerbunt fc ZwaanI (Il98ll ) and iRappaport et al.l (119831 ) was usually adopted. 
However, this law predicts too fast spin-down of lo w-mass MS stars, contradicted with th e 
observation of rapid rotators in young open clusters (ISills et al.ll2000l : lAndronov et al.ll2003l ). 
Obviously a modification of the MB law will have significant infiuence on the period evolu- 
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tion (Ivan der Sluys et al.ll2005bl ). Additionally, there is strong evidence that during LMXB 
evolution the mass transfer is highly non-conservative. Recent measurements of the masses 
of binary and millisecond pulsars indicate that a large fraction of the tra i isferred mass may 
be lo st from the systems rather accreted by the NS (jSassa et al. 2006 : Steeghs &: Jonker 
2007 . and references therein)0. Theoretically , possible ways of ma ss loss have been sug- 



gested, including "evapor ation" of the donor ( Ruderman et al.lll989 ) or " radio-ejection" of 



the transferred material ( Burderi et al.l 12001 



2002 



D'Antona et al. 



20061 ) due to the pul- 



sar radiation/ wind impinging on. In the latter case, the matter is lost from the system at 
the inner Lagrangian (Li) point, carrying away angular momentum and altering the period 
evolution. 

This paper is organized as follows. §2 briefly describes the stellar evolution code, the 
binary models, and the physical assumptions, especially the MB laws and the mass loss 
mechanisms. Then we present the calculated results in §3. Our discussion and conclusions 
are given in §4. 



2. Evolution code and binary model 



2.1. The stellar evolution code 



We use an updated version of the stellar evolution code originally developed by lEggleton 



( ll97ll . I1972I . see also Han et al. 2004, Pols et al. 1995) to calculate the evolutions of binaries 
consisting of an NS (of mass Mi) and an MS secondary (of mass M2). For the secondary 
star we assume a solar chemical composition {X = 0.70, Y = 0.28, and Z = 0.02), the ratio 
of mixing length to pressure scale h eight a = 2.0, and conve ct ive overshooting parameter to 



be .12. The opacity table is from [Rogers fc Iglesiad (119921 ). [Alexander fc FergusonI (119941 ) 
and lHubbard fc Larnpd (119691) . The effective radius of the Roche lobe for the secondary is 
taken from Eggleton ( 19831 ). 



R 



'L,2 



0.49g-2/3 
0.6g-2/3 + ln(l + g-i/3)" 



:i^ 



where q = M2/M1 is the mass ratio, and a is the orbital separation. Mass transfer rate via 
Roche lobe overflow is evaluated as — M2 = RMT ■ max(0, {R2I Rl,2 — l)^)^0yr~^, and we 
adopt RMT = 10^ in the calculations. 



^The massive (~ 1.9Mq) NS discovered in the globular cluster M5 (jFreire et al.ll2008l ) may reflect a bi- 
modal distribution of the initial masses of NSs (rather heavy accreti on during the previo us LMXB evolution), 



as already predicted by hydrodynamical core collapse simulations (jTimmes et al 



199' 



i). 
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2.2. Mass and Angular momentum loss mechanisms 

For LMXBs the timescale of tidal synchronization is much shorter than the characteristic 
evolutionary timescale of the binary, so we can assume that the spin of the secondary star 
and the binary orbital revolution are always synchronized. Assuming rigid body rotation of 
the secondary star and neglecting the spin angular momentum of the neutron star, the total 
angular momentum of the binary system can be expressed as 

J = I2UJ + Jorb 

= l2UJ + G^^^MiM2{Mi + M2)-^^^uj-^^^ (2) 

where I2 is the moment of inertia of the secondary star, u is the angular velocity of the 
binary. 

We consider three kinds of mechanisms of angula r momentum loss. The first is the 



angular momentum loss due to gravitational radiation (ILandau fc Lifshitall975l ) 



dJcR _ 32 G'/^ MlMljM, + M2f'^ 
dt 5 c5 a7/2 ' ^ > 

where c is the light speed. This mechanism is important only in very short period binary 
systems. 

The second angular momentum loss mechanism is for non- conservative mass transfer. 
We assume that a fraction a of the transferred mass is accreted by the NS, and the remaining 
mass is ejected out of the binary as isotropic winds from the NS, carrying away the specific 
angular momentum of the NS, 

—— = -{I- a)M2[Y^) a ^- (4) 

In our numerical calculations we have set a = 0. Alternatively, if the NS is spun up to be 
a millisecond pulsar, its radiation pressure may be strong enough to halt the transferred 
matter at the Li point and quench the accretion. This "radi o ejection" may cause alrn ost all 



the matter from the secondary to be lost from the binary (iBurderi et al.ll200ll . l2002l ). The 
corresponding rate of angular momentum loss is 

—— = -Maa^io; (5) 

where a^i is the distance from the Li point to the center of mass of the binary system. 

The third angular momentum loss mechanism is MB. For a low-mass MS star with deep 
convection zone, stellar winds which are magnetically coupled with the star can deceler- 
ate the stellar spin efficiently, thus carrying away the orbital angular momentum because 
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of tidal synchronization. Th e widely used for n aula f or such MB effect was postulated by 
Verbunt &: ZwaanI (Il98ll ) and iRappaport et al.l (Il983l ) as 



dJ, 



MB 



3.8 X 10~^°M2i?^u;Myncm. 



4, ,3 



dt 



(6) 



However, observations of rapid rotators in youn g open clusters suggest a modification of the 
MB law at high rotation rate (ISills et al.ll200Cll ). 



dJi 



MB 



dt 

d^MB 
dt 



-Ku^ ( -± 

-K0 



& 



ft- 2 / RiX'' (M, 



-0.5 



U > UJ, 



crit ) 



(7) 



where K = 2.7 x 10^^ gcm^s (lAndronov et al.ll2003l ). cjcrit is the critical angular velocity at 
which the angular momentum loss rate reaches a saturated state, given by (IKrishnamurthi et al 
1997h 



'-^crit (t) 



crit0 " 



where t^q^q and are the global turnover timescales for the convective envelope of the Sun 
at its current age and of the secondary star at age t, respectively. They can be calcu- 
l ated by integrating the i nverse local convective velocity over the surface convective envelope 
( Kim fc Demarquelll996l ): 

dr 

(9) 

R, ^ 

where is the radial distance from the center of the star to the bottom of the sur- 
f ace convective envel ope , and v is the local convective velocity from mixing-length theory 
( Bohm- Vitense igssh. Our calc ulation gives rtp,© ~ 28.4 d, slightl y larger than Ttr^.c^, — 13.8 



n 



van der Sluys et a 



d in 

and 

for different values of rtg,© calculated. 



(l2005br). but consiste nt with the results of iKim fc Demarqud (Il996l ) 



Jung fc KimI (120071 ). See lEggletonl (|2006| . page 46) for the discussion of a possible reason 



Following the suggestion of iPodsiadlowski et al.l (120021 ). we also add an ad hoc factor 



exp(-0.02/gconv + 1) if Qconv < 0.02, 

in Eqs. (6) and (7), where gconv is the mass fraction of the surface convective envelop, to 
reduce the MB effect when the convective envelope becomes too small. 
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2.3. Binary models 

To examine tlie influence of mass and angular momentum loss mechanisms on the pe- 
riod evolution, we construct four models with various mass and angular momentum loss 
combinations: (1) model 1 - conservative mass transfer with traditional MB law (Eq. [6]); 
(2) model 2 - conservative mass transfer with saturated MB law (Eq. [7]); (3) model 3 - 
non-conservative mass transfer with mass loss from Li point (Eq. [5]) and saturated MB 
law (Eq. [7]); and (4) model 4 - non-conservative mass transfer with mass loss from the NS 
(Eq. [4]) and saturated MB law (Eq. [7]). In all the four models, the initial NS mass is set 
to be Ml i = 1.4M0, and the initial mass of the secondary M2,i ranges from 0.5 to 2.0Mq. 



3. Numeric Results 

3.1. The bifurcation periods 

Throughout this paper we deflne the bifurcation period Pbif as the initial binary orbital 
period P\ with a zero-age main-sequence (ZAMS) companion star that separates converging 
from diverging systems. We use Pf to denote the fl nal orbital period after the mass transfer. 



Another deflnition of the bifurcation period used by lPodsiadlowski et al.l (120021 ) is the orbital 



period when the Roche lobe overflow just begins, which is expressed as Priof in this paper. 

The results of the bifurcation periods for the four models described in §2.3 are summa- 
rized in Fig. [Hand Tabled! We also draw the minimum initial period Pzams that corresponds 
to a lobe-fllling ZAMS donor star in Fig. [H Several features are noted for the bifurcation 
periods in Fig. [TJ First, the bifurcation periods for all the four models decrease with increas- 
ing initial secondary mass from 0.5Mq to 1.3Mq. Second, in models with saturated MB, 
there exists an upper limit of the initial secondary mass, beyond which no bifurcation period 
exists. This upper limit is in the range ~ 1.2 — 1.3 Mq for model 3, and ~ 1.3 — 1.4 Mq for 
models 2 and 4. Third, comparing the bifurcation periods of model 1 with those of models 
2 — 4 indicates that the MB law plays the most important role in determining the values of 
the bifurcation periods: different MB laws can change the bifurcation periods by as much as 
~ 60%, compared to ~ 14% (see Table [1]) caused by different mass loss mechanisms. 



In Tabled] we also present Priof following iPodsiadlowski et al.l (I20021). In model 1, we get 
Prinf — 18.3 hr for M2^i = IMq, which is in close line with the result of IPodsiadlowski et al. 



(120021 ) (17.7 hr with Z = 0.001, Y = 0.27), where the difference could be explained as the 
difference between the metallicities we used. When we use saturated MB, Priof decreases to 
~ 11 hr. 
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According to the calculated orbital period evolutions, LMXBs can be classified into 
three categories: the diverging systems with Pf ^ P;, the converging systems with Pf <C Pi, 
and the parallel systems with Pf ~ Pi. As an example, we present the calculated results for a 
IAMq + I.OMq binary in model 3, to illustrate the three kinds of evolutionary sequences in 
Fig. [21 The corresponding bifurcation period is found to be 1.25 day, and the initial orbital 
periods are chosen to be Pi = 1.20, 1.25, and 1.40 days, which represent the converging, 
parallel, and diverging systems respectively. 



3.2. Effect of MB and mass loss 



Pylyser fc Savonijd (119881 ) have emphasized the effect of MB on the evolution of LMXBs. 
Comparing the results of models 1 and 2 presented in Table [H we find that the bifurcation 
periods with traditional MB law are smaller (larger) than those with saturated MB law, 
when the initial secondary star mass is less (larger) than O.7M0. 

Our results suggest that mass loss also influences the value of Pbif, though in an less 
important way compared with MB. The bifurcation periods in non-conservative models 3 
and 4 are lower than those i n model 2, in which conserva tive mass transfer has b een assumed 



his r esult is consistent with lvan der Sluys et al.l (l2005bl ). but contradicted with lErgma et al. 



(119981 ) 



It is also interesting to see whether an UCXB can form with saturated MB. For an 
LMXB with an initial orbital period below the bifurcation period, mass transfer is mainly 
driven by the loss of angular momentu m. The orbital period will decrease with the donor 



mass until a minimum period is reached. iPaczvnski fc Sienkiewicz 



i 19811 ) found a minimum 
I showed that minimum 

orbital periods less than 11 min could be reached for binaries with an initial orbital period 
very close to the bifurcati on period if traditional MB is included, but in a time longer than 
the age of the universe. Ivan der Sluys et al.l (l2005bl ) further investigated this "magnetic 
capture" scenario for the formation of UCXBs. Our calculations show that when the initial 
orbital period is close to the bifurcation period, ultra-compact systems (P < 1 h) can indeed 
form with saturated MB, but also in a time longer than the age of the universe. For example, 
for an LMXB with Ms^i = I.3M0 and Pi = 0.46 day in model 4, a final period of Pf = 22 
min can be reached after > 15 Gyr of mass transfer. All the works done by previous authors 
show that a more efficient angular momentum loss mechanism is required to produce UCXBs 
within 13.7 Gyr in this scenario. 
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3.3. Semi-analytical Method 

In this subsection we will try to use a semi-analytical method to understand our numer- 
ical results. First from Eq. (2), we have the following equation 

3 j ^ 3 Ml M2 I M1 + M2 IP , . 

2 J 2^Mi M2' 2M1 + M2 2P' ^ ' 

If we assume a fraction a of the mass lost by the donor is accreted by the NS, i.e., Mi = 
—aM2, we can write the period derivative as 

p = 3j-A(Mi,M2,a)^, (11) 

where 

^ ' ^ Mi(Afi + M2) ^ ^ 

Our analysis is limited to binary evolution with M2 < Mi. In this case it is clearly seen that 
mass transfer increases the orbital period and angular momentum loss decreases the orbital 
period. The bifurcation period is decided by the balance of these two factors. 

Keeping the orbital period unchanged (i.e. P ~ 0), we calculate the maximum mass 
transfer rates for orbital periods from 0.35 day to 0.95 day. This perio d interval covers the 



whole range of the bifurcation periods obtained in this work and in iPodsiadlowski et al. 



(120021 ). We show the calculated mass transfer rates in Fig. [31 and find that they can be 



fitted by an approximate expression as 

- MiiP) ^ 2.73 X lO"*^(P/day)''-^^(^°-")M0yr~^ (13) 

Then we calculate the mean mass transfer rates with constant P and find that they lie 
between M2(P)/2 and M2(P). This means that if we use Eq. f|T3|) to calculate Priof, it will 
deviate no more than ~ 10% from the true value. The mean mass transfer rates here are 
calculated as follows. We fix the binary period P in a constant value in our code, and then 
evolve the donor from its initial mass M2,i to the time when it loses half of its initial mass 
0.5M2,i. This mass transfer process takes a time of T1/2. Then we use 0.5M2,i/Ti/2 as the 
mean M2 for this period P. If we assume that angular momentum loss is dominated by 
saturated MB when P < 10 d, from Eqs. (2) and ([7]) we have 

jL - ^K,,^ ( ^-o.5.4v^^ 2/3 i^i + M^f' 4/3 

Jorb " ^'''^rJ ^mJ ^ G ' M1M2 

~ -6.2 X 10"^^(P/day)"Vr~^- (14) 
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Here we adopt K = 2.7 x 10^^ gcm^s flAndronov et al.ll2003h . u^nt = 2.9 x 10"^ Hz fISills et al. 



20001 ) . Ml = IAMq, M2 = IMq, and replace the radius of the secondary R2 with its Roche 
lobe radius i?L,2 from Eq. (1). Combining Eqs. ([TT1) -(14) with P = we obtain Pj-iof — 12.4 
hr for conservative mass transfer [a = 1), and Priof — 10.7 hr for non-conservative mass 
transfer (a = 0). These values agree well with our numerical results (~ 10.8 — 11.1 hr). If 
the traditional MB law is used, similarly, from Eq. ([6]) we get 



= -3.8 X 10-''M2Rtco'G-'^' ^^\'lf/^' u;'/' 
U ' M1M2 

~ -4.6 X 10"^(P/day)-2/V"\ (15) 



for Ml = 1.4Mq and M2 = IMq. Combining Eqs. ([II])-(13), (15) with P = we get 
Pj-iof — 22.3 hr for conservative mass transfer {a = 1), whic h is about 20% larg e r tha n 



~ 18.3 hr from our numerical calculations and ~ 17.7 hr in iPodsiadlowski et al.l (120021 ). 
The main reason for this difference is that we use the constant value 1.4M0, IMq for Mi, 
M2 in Eq. (12), which should change with time to ~ 2.2Mq, ~ O.2M0. This will decrease 
the coefficient in Eq. (12) and increase the value of Priof by ~ 10 — 20%. For donor mass 
> I.2M0, this will increase the value of PrM by as much as ~ 30%. So it is better to use 
a = instead of a = 1 for donors mass > I.2M0, which could yield more accurate results 
(from our numerical results we find that the deviation of Priof between conservative and 
non-conservative mass transfer is smaller than 10%). 

Using Eq. (11)-(15), we also compute the semi-analytical results of Priof for 0.5 — 2M0 
donors, and compare them with our numerical results of models 1 and 4 in Fig. 4, where the 
semi 4 results are calculated with a = 0. When we calculate the semi 1 results in Fig. 4, 
for the reasons mentioned above and below, we use a = 1 for M2.i < I.IM0 and a = 
for M2,i > I.2M0, where a = 1 should be used. A few points need to be noted for the 
semi-analytical results in Fig. 4. Firstly with the above mentioned equations it is impossible 
to compute the PrM for binaries with M2.i > I.4M0 under conservative mass transfer, since 
both terms in the right side of Eq. (11) are negative when a = 1 and M2,i > 1.4 M0, 
and there will be no solutions for P = 0. We instead adopt a = when M2,i > 1.4M0 
(from our numerical results we find that the deviation of Priof between conservative and non- 
conservative mass transfer is smaller than 10%). Secondly, Eq. (13) is derived only for IM0 
donor star rather donors in the whole mass range (0.5 — 2 Mq), because in the latter case it 
is impossible to find a unified expression of the mass transfer rate like Eq. (13). As seen in 
Fig. 4, the difference between the semi-analytical and numerical results is generally smaller 
than 20% except for donors smaller than O.7M0. The reasons for the big discrepancies when 
^2,i < O.7M0 are discussed in §4. 



- 10 - 



4. Discussion and Conclusions 

Motivated by new ideas about MB and mass loss in LMXB evolution, we liave made 
a systematic investigation on the bifurcation periods in binary models, taking into account 
different MB laws and mass loss mechanisms. We find that the strength of MB is the 
dominant factor in determining the value of bifurcation periods compared with mass loss. 
The stronger MB, the larger the bifurcation periods. This also results in an upper limit for 
the secondary masses beyond which no converging systems exist. 

In our calculations we assume either fully conservative (models 1 and 2) or non-conservative 
(models 3 and 4) mass transfer to constrain the bifurcation period distribution in different 
mass transfer modes. From the expression of A{Mi, M2, a) we always have 

3M9 Mo 
A{M„ M2, 1) = 3 - ^ < A(Mi, M2, 0) = 3 - j^^^^, (16) 

which means that non- conservative mass transfer contributes more to the increase of the 
orbital period than conservative mass transfer. This explains why we generally have a 
lower bifurcation period in non-conservative mass transfer models (models 3 and 4) than 
in conservative mass transfer model (model 2) under the same MB law. The real situation 
may lie between these two extreme cases. For binary systems with donors M2,i ~ 0.5— O.8M0, 
it would take more that 13.7 Gyr before mass transfer begins via Roche lobe overflow. So 
the bifurcation period for these system seems meaningless, unless there exist some unknown 
mechanisms of loss of orbital angular momentum. 

In our semi-analytical analysis in §3.3 we use the condition P ~ to derive the values 
of Priof- This expression seems different from Pf ~ Priof, which is the original definition of 
bifurcation period. We argue here that these two expressions are roughly the same except for 
binaries with M2,i > 1.4M0 under conservative mass transfer {a = 1), the reason of which 
has been given in §3.3. For M2,i < 1.4M0, we find that P/P always scales with P from 
Eq. (11) and (13)-(15). This means that if initially P > (< 0), P/P will become larger 
(smaller) during the evolution, leading to monotonic increase (decrease) of the period, as seen 
in Fig. [21 So for these systems Pf ~ Priof is approximately equivalent with P ~ 0. Several 
rough assumptions in this semi-analytical method contribute to the discrepancies between 
the semi- analytical results and the numerical results in Fig. 4, especially for M2,i < O.TMq. 
First is the use of P ~ as the definition of Priof, which may not work well sometimes. 
Second is the use of Eq. (13), which is most suitable for binaries with M2,j = IMq as pointed 
out in §3.3. Third is the assumption we made that magnetic braking law is the dominated 
mechanism for the angular momentum loss, while the true case is that the MB may not work 
sometimes (for example when the convective envelop is too small). Fourth is that we use 
a constant initial value of Mi i, M2^i for Mi, M2 in Eq. (12) and Eq.(14)-(15), while in the 
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true case Mi, M2 should change with time. This will cause a big problem for a = 1 when 
> IMq, which has been pointed out in §3.3. Fifth reason is the use of coiCTit,© in Eq- (14) 
as the value of Wcrit for all the donors ranging from O.SMq to 1.3M0. At last we conclude that 
(1)P is a fair definition of bifurcation period, and (2)the period evolution during the mass 
transfer phase is in first approximation sufficiently well described by the balance of mass 
transfer and angular momentum loss caused by MB. For these rough assumptions made in 
this semi-analytical method, its results agree with the numerical results only in the leading 
order. 

Our numerical calculations show tha t there is an upper l i mit fo r the donor mass beyond 



which no converging systems will form. iPylyser fc Savonijd (119881 ) found that, in the case 



of Ml i = 4.OM0, there is no converging system existing if M24 > I.TMq, and concluded 
that for any given initial accretor mass there exists a maximum initial secondary mass for 
the formation of converging systems. From our calculations with Mi^i = I.4M0, we find an 
upper limit for the initial secondary mass M2,i between 1.2 and I.4M0 under saturated MB. 
The reason is that for binaries with a MS donor of initial mass > I.4M0, the bifurcation 
period is shorter than the minimum ZAMS period, so that these systems will diverge. For 
traditional MB, this upper limit is > 2M0, beyond the range of donor masses we adopt. 
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by the Natural Science Foundation of China under grant numbers 10573010 and 10221001. 
B.M. thank P. P. Eggleton, Xuefei Chen, and Xiaojie Xu for their help during this work. 
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Table 1. Calculated results of the bifurcation periods for different binary models 



-Pbif (day) Priof (hr) 



M2,i 


Model 1 


Model 2 


Model 3 


Model 4 


Model 1 


Model 2 


Model 3 


Model 4 


O.5M0 


3.20 


4.08 


4.12 


4.01 


31.8 


16.5 


18.2 


13.6 


0.6Mq 


3.02 


3.41 


3.40 


3.35 


25.6 


11.3 


11.1 


10.0 


O.TMq 


2.86 


2.81 


2.79 


2.77 


21.5 


10.0 


9.8 


9.6 


O.8M0 


2.77 


2.41 


2.38 


2.37 


19.0 


10.1 


9.9 


9.8 


0.9Mq 


2.73 


2.14 


2.09 


2.10 


17.9 


10.7 


10.3 


10.4 


1.0 Mq 


2.28 


1.29 


1.25 


1.27 


18.3 


11.4 


10.8 


11.1 


I.IMq 


1.51 


0.63 


0.59 


0.61 


18.9 


11.7 


10.8 


11.2 


1.2Mq 


1.01 


0.55 


0.48 


0.51 


19.1 


11.8 


10.5 


11.0 


I.3M0 


0.86 


0.52 




0.46 


19.2 


12.0 




11.1 


IAMq 


0.90 








20.7 








I.6M0 


0.95 








22.3 








I.SMq 


0.98 








22.8 








2.OM0 


0.92 








21.6 
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Fig. 1. — Bifurcation periods as a function of the secondary mass in an LMXB for the four 
kinds of models described in §2.3. The dotted hne shows the minimum initial period Pzams 
that corresponds to a Roche lobe filling zero-age main-sequence secondary star. 
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Fig. 2.— Period evolutions of an LMXB with M2,i = I.OMq and P; = 1.20, 1.25, 1.40 day in 
model 3. Circle, triangle, and cross mark the onset of the Roche lobe overflow (RLOF), the 
end of the RLOF, and the end of the calculation, respectively. 
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Fig. 3. — Maximum mass transfer rates in an LMXB consisting of IAMq neutron star and a 
IMq secondary at fixed orbital periods from 0.35 d to 0.95 d. Crosses marks the calculated 
data and the solid line represents a logarithmic fit. 
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Fig. 4. — Comparison of the semi-analytical results of Priof with the numerical values for 
model 1 and model 4. The dotted line shows the minimum initial period Pzams that cor- 
responds to a Roche lobe filling zero-age main-sequence secondary star. Here the semi 1 
results are calculated with a = 1 for M2,i < I.IM©, and a = for M2,i > 1.2Mq where 
a—1 should be used. The reasons why we do this are given in the text. 



